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ABSTRACT 

The method of Darboux transformation, which is applied on cnoidal wave so- 
lutions of the sine-Gordon equation, gives solitons moving on a cnoidal wave back- 
ground. Interesting characteristics of the solution, i.e., the velocity of solitons and 
the shift of crests of cnoidal waves along a soliton, are calculated. Solutions are 
classified into three types (Type-IA, Type-IB, Type-2) according to their apparent 
distinct properties. 
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1 Introduction 



The sine-Gordon equation is important for a number of physical and mathematical 
systems where topological solitons are present. In the case of physical applica- 
tions, it is a canonical model describing charge density waves [1], fluxon dynamics 
in Josephson junctions [2], and DNA promoter dynamics [3], for example. In the 
case of mathematical topics, we can mention the embedding of a surface with con- 
stant negative curvature in the 3-dimensional space [4]. In theoretical physics, it 
was studied as a model system of hadrons with the idea that extended particles in 
quantum field theory can be associated with soliton solutions of integrable equations 
[5]. The quantum sine-Gordon model is also developed and the exact S-matrix of 
the solitons was derived by using the bootstrap methods [6]. 

The important property of the sine-Gordon equation in these studies is the exis- 
tence of soliton solutions. The simple functional form of N-soliton has led to a wide 
range of scientific applications. Considering the potential applicability of the soliton 
solutions, it is desired to have a more generalized form of soliton solutions, having 
more parameters and characteristics. One possible scheme in this direction is con- 
structing soliton solutions lying on a background field. It would give more freedom 
in controlling and characterizing solitons by using proper backgrounds. Indeed, the 
sine-Gordon equation, as an integrable equation, has solutions of superposed states 
"soliton + cnoidal wave". Recently, this type of solutions have attracted new in- 
terest in nonlinear optics, described by the nonlinear Schrodinger equation (NLSE), 
due to the possibility of dynamically reconfigurable photonic structures [7, 8, 9]. 

In fact, the "soliton + cnoidal wave" solution was obtained from the general 
quasi-periodic solutions of N-phase theta functions [10, 11], by taking degenerate 
limit of the 2-phase solution [12, 13]. To apply solutions from the N-phase theta 
functions to real physical problems, the authors in [12, 13] solve the so-called "ef- 
fectivization" problem, which is related to extracting physical solutions by taking 
proper initial conditions [12, 13, 14, 15]. However we find that above solutions, 
though explicit, are inconvenient to analyze their properties, as the parameters in 
the solutions are related in a complicated form. In fact, the mathematical and com- 
putational complexities of computing N-phase theta functions have hindered the use 
of these solutions in applications except some clever calculations like in [16, 17]. 

In this paper, we employ a simple, but powerful soliton finding technique [8, 9] 
based on the Darboux transformation (DT) to find solutions of "soliton + cnoidal 
wave". Compared to the results from the N-phase theta functions, the form of 
solution, especially the parameters, are simple and convenient to analyze various 
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properties of the solutions. Explicit real solutions are described by a parameter u 
without any further constraints and it is directly related to the characteristics of 
solutions. This DT technique gives a clear classification scheme of the solutions 
(Type-lA, Type-IB, Type-2) according to their apparent distinct properties like 
the modulating form of cnoidal waves and the velocities of solitons. The velocity of 
solitons are expressed in a form of relativistic addition, which is a unique feature of 
the "soliton + cnoidal wave" solution of the relativistic sine-Gordon equation and 
is essentially different from that of the nonrelativistic NLSE and KdV equation. As 
the sine-Gordon theory has various physical applications, our solution could yield a 
new interest on controlling the soliton using background fields. 



2 sine-Gordon equation 

2.1 Cnoidal wave solutions 

The sine-Gordon equation is 

dd(p = 2(3 sin 20, (1) 

where d — -§^,d — Jj, and z = t — x and z = t + x represent the lightcone 
coordinates, while x, t are ordinary space-time coordinates. It has two types of 
cnoidal wave solutions; 

$\z, z) = sin- 1 (sn( X , k)), 8^ = \\[^Mx, k), Type - 1 
^\z,z) = sm-\ksn(k X ,k)), 00« = 2kMcn(k X , k), Type-2 (2) 




where 

X = 

and sn, dn are the standard Jacobi elliptic functions. V is related to the velocity 
v c of the cnoidal wave as V = (1 + v c )/(l — v c ) and k G (0, 1) is the modulus of 
the Jacobi function. Here we consider the case of V > 0,\v c \ < 1 and (3 > 0. The 
V < 0,\v c \ > 1 case, corresponding to the rescaled inverted pendulum [17], can be 
obtained by considering the x <-> t, (3 <-> — (3 symmetry of the sine-Gordon theory. 
As far as elliptic functions are involved we employ terminology and notation of Ref. 
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[18] without further explanations. 0^ of the Type-1 solution is a monotonically 
increasing or decreasing function of x, while (p^ of the Type-2 is an oscillating 
function of X- 



2.2 Lax pair 

To obtain a superposed "soliton + cnoidal wave" solution using the DT method 
[19, 20, 21], we need to find a solution of the following linear equations associated 
to the sine-Gordon equation (Lax pair), 

(d + ip\) Sl + i{d<j>f)s 2 = 0, (d - i(3\)s 2 + i{d<f>f)sx = 0, 

(9 _ * cos(2#)) Sl + 1 sin(20«) S2 = 0, 

(5+1 cos(20«))s 2 - i sin(20«) Sl = 0, (4) 

where A is an pure imaginary number and % — 1,2 represents two types of cnoidal 
solutions in Eq. (2). We note that the compatibility between the equations (4) 
requires <f$ satisfies the Eq. (1). 



2.3 Sym's solution 

The solution of the linear equations (4) for the cnoidal wave ((j>^\i = 1,2) can be 
obtained by a slight modification of Sym's solution [22]. Sym first introduced this 
type of solution in a context of vortex motion in hydrodynamics. It was then applied 
to NLSE- related problems in Ref. [23, 8, 9]. A more detailed proof of Sym's solution 
(in a slightly different notation) is given in the Appendix of Ref. [9]. The Sym's 
solution in the case of sine-Gordon equation for the Type-1 cnoidal solution, (f>^\ is 
(to be called as the Type-lA solution) 

- = {^e, (x - u) + ft ^ -<y + U| kf »e,( x + » )} /e.(x), 

where a, 6 are arbitrary real parameters. Here, 

M= /^ cnKA:) 4 + (A: 2 -l)snKA:) 4 5 | /9'» | dn(-u, fc)(l — 2sn(tt, fc) 



/csn(w, fc) cn(w, fc)dn(u, k) \Q s (u) 2sn(w, k) cn(w, k) 
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and 

©.(«) = OoQ) = 1 + 2 E(-) V 2 cos(^), (7) 

with g = exp(—nK'/K) and Jacobi theta function 9 . The parameter -u is related 
to A as 

A = i dn(M, k) 

\JV ft k cn(w, k) sn(w, k) 

3 Type- 1 A soliton solution on a dn-type back- 
ground 

3.1 soliton solution 

A new solution C _ S , describing a soliton moving on a cnoidal wave <f>^\ can be 
constructed using the DT method [19, 20, 21] as following; 



s 1 s 2 

sin(20 c _ s ) = 4 r i ;V" 1 '7 cos(20W) - ^ .' ^ sin(20( 1 )). (9) 



-(2^)) - gi±|gM±j s ,,,, ,U ^ ■ 
(sf-s 2 .) 2 (s?-s|) 2 



Using the fact that Sj,i = 1, 2 in Eq. (5) satisfy the associated linear equations (4), 
it can be proved that <fi c -s in Eq. (9) satisfies the sine-Gordon equation (1). It can 
be seen that the reality of C _ S requires u, as well as a, b, are real in the Type-lA 
solution. 

Figure 1(a) plots <90 c _ s as a function of time t and space x. It is obtained by 
using Eqs. (9), (2), (5) and (6). It shows a characteristic soliton of sine-Gordon 
equation lying on a cnoidal wave background. The parameters used are V = 1, k — 
0.8, u = 0.7, (5 — 1, a — b — 1. Figure 1(b) shows a time sliced view of the soliton in 
Figure 1(a) at t — 0. These figures are drawn using MATHEMATICA, which is also 
used to check that the solution in Eq. (9) indeed satisfies the equation of motion 

(!)• 

3.2 Properties of the solution 
3.2.1 Shift of cnoidal crests 

In figure 1, we can see that C _ S becomes a cnoidal wave when we move away from 
the soliton. In fact, as M — > oo, s 2 — > — isn(u, k)dn(x— u, k)si/cn(u, k) and <90 c _ s — > 
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Figure 1: <90 c _ s showing a sine-Gordon soliton residing on a cnoidal wave 
background (Type-lA solution); (a) 3-dimensional plot, (b) time-sliced view 
at t — 0. The parameters are V — 1, k — 0.8, w = 0.7, f3 — \,a — b— \. 



-oo, 



-2yj0/Vdn(x-2u, k)/k, which is a cnoidal wave. Similarly, in the region M - 

s 2 — > — icn(/u, fc)si/(dn(% + /c)sn(u, k)) and (90 c __ s — > — 2^J f3/Vdn(x + 2m, 
This calculation shows that crests of the cnoidal wave are shifted by Au across the 
soliton, which can be seen in Fig. 1. Explicitly, the shift of crests for parameters of 
Fig. 1 is 4m = 2.8 in x and ku = 0.56 (V = 1, v = 0, (3 = 1) in x, use Eq. (3). 



3.2.2 velocity of the soliton 

The soliton in Fig. 1 moves along a line described by M = 0, which can be casted 
in a form = z + Az with 

A _ y | y cn(u,ky+(e-l)sn(u,ky ,( 2 Q' s (u) | dn(M,A:)(l-2sn(M,A:) 2 )^ 



sn(w, k) cn(w, k)dn(u, k) 



e s (u) 



sn(u, k) cn(w, k) 



(10) 

Using the space-time coordinates t, x instead of the lightcone coordinates z, z, it is 
written as = z + Az = 2(x — vt)/(l — v) where v — (A + 1)/(A — 1). It shows v is 
the velocity of the soliton in the ordinary space-time coordinates. We can express 
the soliton velocity v in terms of the v c (V = (1 + v c )/(l — v c )) of the cnoidal wave 
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Figure 2: Intrinsic velocity v versus u, Type-lA sine-Gordon soliton for 
k = 0.5 (solid line) and k = 0.9 (dashed line). 



as v — (v c + v o)/(l + v c Vq) where 

A; 2 sn 2 w cn 2 w + dn 2 -u (1 — 2sn 2 w) 

= v(v c = 0) 



— 4snu cn-u dnu Q' s (u)/Q s (u) + k 2 sn 2 u cn 2 u — dn 2 u (1 — 2sn%)' 

(11) 

where sn-u = sn(u,k), etc. It shows that vq is the intrinsic velocity of the soliton 
observed in the moving frame of the cnoidal wave, which depends on the parameters 
u and k. The soliton velocity v is given by the relativistic addition of two velocities, 
v c (velocity of the cnoidal wave) and v (intrinsic velocity of the soliton), where v c 
does not depend on k or u. Figure 2 shows the intrinsic velocity v in u for a soliton 
on cnoidal wave backgrounds of k — 0.5 (solid line), and k = 0.9 (dotted line). 
Note that vq is a periodic function in u with the periodicity K. (K(k = 0.5) = 
1.69, K(k = 0.9) = 2.28) It shows that the intrinsic velocity vq — — 1 at u — nK for 
integer n. Using Eq. (11) and the relation 

e ^~ , ^=W, + ^ u ' (i2) 

we can find that the velocity v = —k'K/E at u = K/2. Especially, this value is 
-0.995 for k = 0.5, -0.848 for k = 0.9, -1 for k = 0, and for k = 1. Thus a 
soliton can be stopped only when it resides on a cnoidal wave with k — 1. On a 
k = background, the soliton moves with a velocity t>o = — 1. Note that there is 
a symmetry, t — > — t, in the "soliton + cnoidal wave" solution in Eq. (9). This 
symmetry gives a solution having a soliton velocity —Vq. 
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4 Type- IB soliton solution on a dn-type back- 
ground 



4.1 soliton solution 

The reality requirement of C _ S was satisfied by taking real u in the case of Type 
-1A solution. There is another solution for the reality requirement, which is by 
taking u = w + iK', with real w (to be called as the Type-IB solution). In the 
Type-IB case, we obtain following expressions by substituting u = w + iK' in the 
corresponding expressions of the Type- 1 A solution; 

o / / — \ 2 [P A f n . Al [P~sn(w, k) cn(w, k) S 

where 

acn(x — w , k)e M Q s (x — w — iK')/sn(x — w,k) — bdn(w, k)e~ M Q s (x + w + iK') 
bcn(x + w, k)e~ M Q s (x + w + iK')/sn(x + w,k) + adn(w, k)e M Q s (x — w — iK') ' 

(14) 

and 

M = _ ftp dn( W , kf + k' - 1 _ z ( e' s {w + iK') + cn{w,k){k\n{w,kf-2) \ ^ 
* A; 3 sn(w, /c) cn(w, /c) dn(w, A;) \O s (w + ifT') 2dn(w, A;)sn(w, k) ) 

' (15) 



4.2 Properties of the solution 

Figure 3(a) plots <90 c _ s of Type-IB as a function of time t and space x. The pa- 
rameters used are those of Fig.l, i.e., V — 1, k — 0.8, u = 0.7, (3 — 1, a — b — 1. The 
distinct difference between the Type-lA and Type-IB solution can be seen from a 
time sliced view of the soliton at t — 0, see Fig. 3(b) and compare it with Fig. 1(b). 
The intrinsic velocity v of the Type-IB soliton can be obtained from that of the 
Type-lA soliton by substituting u — * w+iK', which gives v (u = w+iK') = —v (w). 
Thus the characteristics of the intrinsic velocity vq of the Type-IB soliton is essen- 
tially the same as those of the Type-lA soliton. The shift of crests of the cnoidal 
wave can be similarly calculated as in the Type-lA solution, which is Aw. 
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Figure 3: <90 c _ s showing a sine-Gordon soliton residing on a cnoidal wave 
background (Type-IB solution); (a) 3-dimensional plot, (b) time-sliced view 
at t — 0. The parameters are V — 1, k — 0.8, u = 0.7, f3 — \,a — b— \. 



5 Type-2 soliton solution on a cn-type background 
5.1 soliton solution 

The "soliton + cnoidal wave" solution with the ^ background can be obtained 
by substituting k — > l/k,u — > ku,x ~^ k 2 x in the Type-lA solution. We call it 
the Type-2 solution. Note that dn(%, k) — > d\i(k 2 x, 1/fc) = cn(/cx, fc) such that 
d(j)^(z,z) — > d(j)@)(z,z). In this case 

90 c _ s (z, *) = 2kMcn(k X , k) + 4</f C ; (M ; fc) rr^r , (16) 
V V V v sn(u, fc) dn(w, k) S 2 + 1 



(17) 



where 

_ aA;sn(w, A;)cn(A;x - u, k)e M Q t (X2 - u) + bdn(u, k)e~ M Q t {x2 + u) 
bksn(u, k)cn(kx + u, k)e^ M Q t (x2 + u) + adn(-u, k)e M Q t (x2 — u) ' 
with 

M _ r^ l -2k 2 sn{u,k) 2 + k 2 sn{u,k) A ^ | U &M | kcn(u, k)(l — 2k 2 sn(u, k) 2 ) \ 
* sn(w, fc) cn(w, k) dn(u, fc) \ ©t(«) 2 dn(w, fc) sn(w, fc) / 

(18) 

and 

e <<") - e ° {m^)) = 1 + 2 ^ ( - Tq " cos( ^ ) - (19) 
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Figure 4: <90 c _ s showing a sine-Gordon soliton residing on a cnoidal wave 
background (Type-2 solution); (a) 3-dimensional plot, (b) time-sliced view 
at t — 0. The parameters are V = l,k = 0.8, w = 0.7, f3 — \,a — b— \. 



with g = exp[-7r^7(K - iK')\. 

5.2 Properties of the solution 

Figure 4(a) plots <90 c _ s of the Type-2 soliton as a function of time t and space x. The 
parameters used are those of Fig.l, i.e., V = 1, k = 0.8, u = 0.7,(3 — l,a — b — 1. 
The distinct difference between the Type-2 and Type-1 solutions can be seen from 
a time sliced view of the soliton at t — 0, see Fig. 4(b) and compare it with Fig. 
1(b) and Fig. 2(b). The intrinsic velocity vq of the Type-2 soliton can be obtained 
from that of Type-lA soliton by substituting k — > 1/k, u — > fcu, which gives 



Vq = 



dn 2 u sn 2 u + cn 2 w (1 — 2A; 2 sn 2 -u) 



sn-u cnu dn-u(— AQ' w (u) /Q w (u) + 2nu/(KK')) + dn 2 w sn 2 -u — cn 2 -u (1 — 2k 2 sn 2 u) 
where we use the identity 

©{(«) 



u 



e t («) 



— 7T 



+ 



2{K-iK')K> e w (u) 



(20) 
(21) 



with 



e.(«) = feH^) = 2^ £ ^ cos( -'^+ 1 )- ), 



2X' 



(22) 
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Figure 5: Intrinsic velocity v versus u of the Type-2 sine-Gordon soliton for 
k = 0.5 (solid line) and k = 0.9 (dashed line). . 

with q = exp(—nK/K'). The characteristics of the intrinsic velocity v of Type-2 
soliton is shown in Fig. 5. Special values are vq = — 1 at u = 2nK, n=integer and 
v o = 1 at u = (2n+l)K. At u — (n+l/2)K, vq = ±oo. It shows that type-2 solitons 
move with the velocity \v\ > 1. We can obtain soliton solutions having \v\ < 1, by 
taking the Backrund transformation (0 C _ S , x, t) <-> (0 c „ s + 7r, t, x). The shift of crests 
of the cnoidal wave can be similarly calculated as in the Type- 1 A solution, which 
gives 4m in the Type-2 solution. The substitution u = w + %K' (used to obtain the 
Type-IB solution from the Type-lA solution) does not give new solution in the case 
of Type-2 solution. In fact, there exists a symmetry u <-> —u + K + %K' in the 
Type-2 solution. 

In the k — > limit of the Type-2 solution, 4>^ — > 0. Using the fact that 
K(k = 0) = tt/2, K'(0) = co, we obtain S = - exp(-2M + In b/a) where 

1777, ( sin u cos u \ cos u , . . 

v Vcosm sinu/ 2sm-u 

Especially at w = if/2 = tc/A with a = 6 (f3 = V = 1), we obtain 

C _ S = -2tan _1 exp(4t), 90 c _ s = -2sech4t, sin20 c _ s = 2sech4t tanh4t. (24) 

Considering the x <-> t, (3 <-> — /3 symmetry of the sine-Gordon theory, the solution 
in Eq. (24) becomes the well-known sech-type soliton (without a background) of 
P — — 1 theory. 
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6 Discussion 



In this paper, we introduce "soliton+cnoidal wave" solutions of the sine-Gordon 
equation. It was obtained using the DT and the Sym's solution of the associated 
linear problem on a cnoidal wave background. The solutions are expressed in terms 
of the Jacobi elliptic functions and can be easily manipulated to obtain interesting 
physical characteristics. In fact, we use the symbolic package MATHEMATICA to 
check various formulae of the present paper including the solutions itself, as well 
as to obtain the figures. We calculate the velocity of a soliton on a cnoidal wave 
and the shift of crests of the cnoidal wave. The velocity of the soliton is given by 
the relativistic addition of two velocities, the velocity of the cnoidal wave and the 
intrinsic velocity of the soliton observed in the moving frame of the cnoidal wave. 
The intrinsic velocity is determined by the property of the cnoidal wave (dn or cn 
wave and the modulus k) and the DT parameter u. It would be interesting to make 
an (numerical) analysis about the nature of the solitons, especially when it lies on 
a finite-width cnoidal wave background. 

The stability analysis of these solutions is remained for future study. There 
already appeared that linearized instabilities of N-phase solutions can be labeled 
in terms of spectral data [17]. It implies that there exist regions in the parameter 
space of u, k and L (finite length of cnoidal waves), where the solutions are stable. 
Explicit relation between the spectral data and DT parameter u is not understood 
yet. 
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